library("Rmpi")
library("doMPI")

cl <- startMPIcluster(count=2)
registerDoMPI(cl)
foreach(i=1:100) %dopar% sqrt(i)
closeCluster(cl)
